A Line of Critical Points in 2 + 1 Dimensions: Quantum Critical Loop Gases and Non-Abelian 

Gauge Theory 
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We (1) construct a one-parameter family of lattice models of interacting spins; (2) obtain their exact ground 
states; (3) derive a statistical-mechanical analogy which relates their ground states to O(n) loop gases; (4) show 
that the models are critical for d < where d parametrizes the models; (5) note that for the special values 
d — 2 cos(7r/(fc + 2)), they are related to doubled level-fc SU(2) Chern-Simons theory; (6) conjecture that they 
are in the universality class of a non-relativistic SU(2) gauge theory; and (7) show that its one-loop /3-function 
vanishes for all values of the coupling constant, implying that it is also on a critical line. 
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Family of Microscopic Lattice Models. Consider the fol- 
lowing models of s = 1/2 spins on the links of a honeycomb 
lattice, inspired by a model of Kitaev |l|]: 
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where A/"(u) is the set of 3 links neighboring vertex v, and 



F° 
p 

F 1 

F 2 
1 p 

F 3 
1 v 



a 1 o~2 cr 3 <r 4 a 5 a 6 
a i a 2 (T z (J i' J ^' J Q + cyclic perm. 



0", <7 <7, &4 0"s cr. 



cyclic perm. 
°i cr 2 _(T 3" (J 4~ cr r cr 6^ + cyclic perm. 



(2) 



if 1, 2, . . . , 6 label the six edges of plaquette p. This Hamilto- 
nian is a sum of projection operators with positive coefficients 
and, therefore, is positive-definite. The eigenvalues of the first 
term in Q are 2, 0, corresponding to whether there is an even 
or odd number of cr 2 = — 1 spins neighboring this vertex. The 
zero eigenvalue corresponds to the latter case. When eigen- 
value zero is obtained at every vertex, the a z = 1 links form 
connected loops. Hence, the zero-energy subspace of the first 
term is spanned by all configurations of multi-loops. Loops 
on the honeycomb lattice cannot cross. 

The term on the second line is a projection operator which 
annihilates a state \ $>) if the amplitude for all of the spins on a 
given plaquette to be up is a factor of d times the amplitude for 
them to all be down, i.e. if the amplitude for a configuration 
with a small loop encircling a single plaquette is a factor of 
d times the amplitude for an otherwise identical configuration 
without the small loop. The other three lines of the Hamil- 
tonian vanish on a state | if it accords the same value to a 
configuration if a loop is deformed to enclose an additional 



plaquette. It is useful to think of these multi-loops as continu- 
ous curves, in which case, the ground state is invariant under 
smooth deformation of a curve and it loses a factor of d if a 
small, contractible curve is erased. Together, these conditions 
have been dubbed '<i-isotopy' J3l. 

Exact Ground State Wavefunctions. The terms on the final 
four lines of Q do not commute with each other, but they 
are compatible in the sense that they all annihilate the ground 
state, which is a superposition of all configurations a of multi- 
loops weighted by a factor of d to the number of loops n a in 
each configuration: 
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On an L x L torus, the ground state degeneracy is ~ L 2 be- 
cause the Hamiltonian does not mix states \a) with different 
winding numbers. The different ground states are given by 
(0 but with the sum over a restricted to a single topological 
class. More generally, the ground state on any genus g > 1 
surface is infinitely degenerate in the thermodynamic limit. 
At d = 1, there is an additional condition called the Jones- 
Wenzl projector which also annihilates the ground state (TJ} on 
a topologically-trivial manifold but mixes different winding 
number sectors on higher-genus surfaces. Hence, at d = 1, 
there are two Hamiltonians which have the same ground state 
on the sphere but one of them, given in eq. ^ is infinitely de- 
generate on the torus while the other (with the Jones-Wenzl 
projector added) has finite degeneracy 01 • The second leads 
to a topological phase with an energy gap 1 1], while the first, 
as we will see later, is gapless and critical. At intermediate 
length scales, the two systems are the same and their physics 
is determined by d-isotopy; it is only at longer scales (where, 
for instance, the topology of the manifold becomes apparent) 
that, in the latter Hamiltonian, the critical behavior crosses 
over to that of the stable phase. 

While the ground state can be obtained exactly, excited 
states cannot because the different projection operators on the 
final four lines of do not commute with each other. We 
will obtain some information about excited states using a vari- 
ational ansatz but, in order to do this, we need to learn a little 
more about the structure inherent in the ground state. 
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Mapping of the Ground-State to a Statistical-Mechanics 
Problem. Many properties of the ground-state wavefunction 
can be obtained by observing that the norm of the ground 
state is equal to the partition function of a classical 0(n) loop 
model with n = d 2 : 

<*o|*o) = E^"" =Zo(n)(x = n) (4) 

a. 

For integer n, this model can be defined by the partition func- 
tion 

Z (n){x) = f Y[dSi J] (1 + xSi ■ 4) (5) 

where Si lies on the unit (n — l)-sphere. The Hamiltonian 
— j3H = j) m (l + x $i • Sj) has been chosen so that its 
high-temperature (small x) series expansion takes the form 

3o W (*)=£(jir ^ ( 6 ) 

a 

where £ a is the total length of the loops in the configuration a. 
The expression l|6} is well-defined for arbitrary n and x, so we 
will take it as the definition of the 0(n) loop model y|. For 
n < 2 and x > x c = n\/2 + \/2 — n, this model is in its low- 
temperature phase, which is critical. Spin-spin correlation 
functions have power-law decay 0|, (S(r) S(0)) ~ r~' n , 
where i]k — \k 2 g — |(1 — g) 2 and < g < 1 is given 
by n = —2 cos(7rg). In this regime, the loops meander about 
the system, as described by a family of exponents such as r\. 
(These exponents are obtained for any oo > x > x c , and 
our results apply to the corresponding family of wavefunc- 
tions within this universality class. Precisely the same expo- 
nents are also obtained in the critical g-state Potts model with 
q = n 2 on the square lattice 01, where loops surround 
clusters with equal Potts spin. This can be exploited in con- 
structing other lattice models in the same universality class 

0.) 

However, equal-time spin-spin correlation functions such 
as (erf cr| ) in the original quantum model Q are short-ranged 
in space. Such correlation functions are related to the prob- 
ability that a loop passes through i and a loop which may or 
may not be distinct passes through j. Such correlation func- 
tions vanish in the 0(n) loop models (or the related g-state 
Potts models). These models have a Coulomb gas representa- 
tion as Gaussian height models with background charges |3] 
whose presence causes correlators of neutral operators, such 
as gradients of the height (to which the local loop density cor- 
responds) to vanish. Algebraic decay is possible for correla- 
tion functions of operators which are charged in the Coulomb 
gas picture, but these are non-local in terms of the spins of 
since they measure, for instance, the probability that two spins 
of and cj lie on the same loop. (At the two points d = 1, \/2, 
this can also be seen from the fact that the ground state on the 
sphere is the same - and, therefore, has the same equal-time 
correlation functions - as that of a gapped Hamiltonian |(J 



which is a sum of local commuting operators and, therefore, 
has correlation length zero.) Thus, the ground state wavefunc- 
tion of Q has an underlying power-law long-ranged structure 
which is apparent in its loop representation, but it is not mani- 
fested in the correlation functions of local operators of. As we 
will see momentarily, this long-range structure leads to gap- 
less excitations for the Hamiltonian Q and, therefore, long- 
ranged correlations in time in spite of the lack of long-ranged 
correlations in space. We call such a state of matter a topolog- 
ical critical or quasi-topological phase. 

Low-Energy Excitations. In spite of the short-ranged na- 
ture of equal-time spin-spin correlation functions and the ab- 
sence of any conservation laws for the Hamiltonian Q, we 
can construct a variational argument that this Hamiltonian is 
gapless using the criticality of non-local correlation functions. 
The key observation is that the configuration space can be 
divided into two regions X, Y - those configurations which 
have 'large' loops (X) and those which don't (Y). Let us 
define Targe' to mean having a linear extent (the greatest dis- 
tance between two points on the loop) which is larger than uL, 
where L is the linear size of the system and < u < 1 will be 
choosen in a moment. If we consider a large but finite-sized 
system, there is a finite probability p(u) that there will be a 
large loop so long as 1 < d < \/2. Since the loop model is 
critical, p can depend only on the shape of the system and u, 
not on its size; this is the key input following from the critical- 
ity of the 0(n) loop model. We choose u so thatp(w) = 1/2 
and consider a trial wavefunction: 

1*1) = -j^= ( i«) - E rf,lQ i«> ) ^ 

The prefactor 1 / W 'Zo(n) normalizes the wavefunction. The 
condition p(u) = 1/2 guarantees that (^ol^i) = 0. To com- 
pute the energy of this wavefunction, we note that the opera- 
tors Fp' 2 3 can change the spatial extent of a loop by at most 
one plaquette, i.e. by less than 2a where a is the lattice spac- 
ing. Hence, the Hamiltonian -ffd-iso nas the special property 
that (a'\H\a) — for a £ X and a' £ Y unless a' S dY 
(this is necessary but not sufficient), where dY is the set of 
loop configurations whose largest loop has linear extent R sat- 
isfying uL — 2a < R < uL. If the matrix element is non-zero, 
it takes the value —1. These considerations, along with the 
fact that ffd-isol^o) = 0, tell us that 

(*i|ff d -u»|*i) = 

-T^E E^"" (*'\Hd-iso\a) 

<-A- E d2< = 4(p(u-2a/L)-p(tO) 

Z 0(n) a , edY 

«4p»^ (8) 

Hence, the energy gap vanishes in the L — ► oo limit for 
1 < d < \/2. The bound (|8j can be tightened by noting that 
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this eigenvalue problem is analogous to that for a stretched 
string; by optimizing the trial wavefunction, a more care- 
ful treatment shows that the gap vanishes with system size 
as 1/L 2 . Details will be given in ref. 7, where we will 
also present a more general theorem which applies to Hilbert 
spaces which decompose into two subspaces X, Y and Hamil- 
tonians which, like iJd-iso, do not directly connect states in 
these two subspaces, but connect them only through a long se- 
ries of repeated applications of it. At some intermediate step 
in this passage from X to Y, the system must pass through a 
bottleneck; the configuration space has the 'dumbbell' form 
depicted in figure ^ In our case, the bottleneck corresponds 
to those configurations in which the longest loop has length 
uL. From this theorem and the analogy drawn in figure^ we 
see that the eigenvalue problem for our Hamiltonian is analo- 
gous to the eigenvalue problem for the Laplacian operator on 
a manifold shaped as in figure U The square lattice quantum 
dimer model at its RK point ||8|] is another example of such a 
system, with X and Y corresponding, respectively, to config- 
urations with large (X) and large negative (Y) total areas (or 
integrated heights, in their associated height models) enclosed 
by the curves in their transition graphs. 

From the size-dependence of the gap, we deduce that the 
low-energy excitations of i/d-iso have dispersion u oc fc 2 . 
However, these gapless modes are not Goldstone bosons, 
since there is no continuous symmetry of Q. They are crit- 
ical modes which are unstable to the addition of terms to 
our Hamiltonian which allow surgeries which cut and rejoin 
curves, such as the Jones-Wenzl projectors |2]. Such terms 
would directly connect the left and right sides of figure \l\so 
that the configuration space of the system would have the 
shape of a ball instead of a dumbbell and, hence, would be 
gapped. 



config. without 
large loops 




config. with 
large loops 



FIG. 1: The configuration space of the system has the shape of a 
dumbbell with two regions X, Y which are not directly connected by 
the Hamiltonian. Consequently, a trial wavefunction which has value 
-$o on y and $o on X has vanishing energy in the thermodynamic 
limit. 

Conjectured Effective Field Theory. We do not have a di- 
rect derivation of the low-energy effective field theory for this 
line of critical points, but we can motivate a conjecture with 
the following reasoning. Quantum loop gases correspond nat- 
urally to Wilson loops in gauge theories. The gauge group 
should be SU(2) because its fundamental representation is 
pseudo-real which implies - as in the Rumer- Teller- Weyl the- 
orem - that the loops are undirected. Quantum gases of undi- 
rected loops arise in the context of doubled SU(2) Chern- 
Simons topological field theories |2], with which our critical 
points are clearly intimately related. The effective field theory 



should have extensive degeneracy on the torus corresponding 
to the different distinct winding numbers of the loop gas. Fi- 
nally, the effective field theory should have quadratic disper- 
sion u) oc fc 2 , rather than the linear dispersion associated with 
relativistic gauge theories. The only natural choice for such a 
theory is a variant of SU(2) Yang-Mills theory with the coeffi- 
cient of the usual electric field term tr(E 2 ) set to zero (the ab- 
sence of such a term implies the existence of infinitely-many 
degenerate ground states on the torus, corresponding to dif- 
ferent uniform electric fields) so that the leading electric field 
term has two derivatives |9]: 

S[A a ,A^,Ef] = ±J d 2 xdT^Efd T A? + A a D t Ef 
- \ {DjEff + l -B a B a + \x{E?EZf + \ 2 {EfEj) 



(9) 



where D~E? 



djEf + f abc A b Ef and the magnetic field is 
given by B a = d Y A a 2 - d 2 A\ + f abc A\A c 2 . The SU(2) index 
a takes the values 1,2,3, while the spatial indices i, j = 1,2. 
jabc are tne structure constants of SU(2). Af , Ef are an SU(2) 
gauge field and its canonically conjugate electric field. The 
time component of the gauge field, Aq, is a Lagrange multi- 
plier which enforces Gauss' law, DiEf = 0. We have written 
the action in first-order phase space form. In principle, Ef 
can be integrated out so that we will have an action dependent 
on A1 alone, but this is cumbersome because Ef and d T A? 
are not linearly-related, unlike in ordinary Yang-Mills theory. 
The theory is invariant under the usual gauge transformation 
A l + ^Q a + f abc A b a c , Ef -> Ef + f abc A b a c . Two 
quartic terms, with couplings Ai, \ 2 have been included be- 
cause they are marginal at tree-level. All other higher-order 
terms are irrelevant or forbidden by symmetry and have been 
dropped. 

RG Analysis of Non-Relativistic Non-Abelian Gauge The- 
ory. This theory is marginal at tree-level, as a result of its 
uj oc k 2 dispersion, which raises the concern that it might be 
massive as a result of quantum fluctuations, just as AD Yang- 
Mills is. To address this, we compute the /3-function for this 
theory to one-loop using the background field method and di- 
mensional regularization . Remarkably, it vanishes: 



dg 
di 



^7=0 



(10) 



The cancellation appears to result from the similarity between 
the way in which the gauge field Af and its conjugate momen- 
tum Ef enter the first four terms of the action (|9jl, which gives 
us hope that it survives to all orders. (For details, see ref. 0.) 
For the same reason, the relative scaling of space and time, 
which can flow in principle, is unrenormalized at the one-loop 
level, so that the dynamical exponent remains z = 2. Since 
the theory (|9jl is on a critical line for A1.2 = 0, it is a viable 
candidate theory for the universality class of Q. However, 
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the one-loop RG equations for Ai,2 generically run away: 

^ = 9 — (24Ai + 20A 2 - 22A? - 28AiA 2 - 
dt 16 2 

^ = -fg ( 8A i + m 2 + 4A? + 16AxA 2 + 14A|) (11) 

so that (|9jl is not critical for small g, where a one-loop cal- 
culation can be trusted. This suggests the following picture, 
which echoes our earlier analysis of Q. The theory is gapped 
for small g as a result of the runaway flow of Ai j2 ; this cor- 
responds to the regime d > \/2 where Q is presumably 
gapped. On the other hand, we expect that for g sufficiently 
large, Ai, 2 will be irrelevant (this needs to be checked by a 
higher-loop calculation of the RG equations for Ai j2 ), so that 
l|9} is critical, corresponding to the critical regime d < \/2. 
A tr(E 2 ) = EfEf term is relevant at g = 0, so its coeffi- 
cient must be tuned to zero. It is an important open question 
how relevant this term is at the large g which we expect to 
correspond to d < \/2. 

Correlation Functions. At present, we do not know how 
to directly compute correlation functions in the critical gauge 
theory l|9). However, if the correspondence with the <i-isotopy 
critical theories is correct, we can deduce some non-trivial 
equal-time correlation functions involving Wilson loop oper- 
ators WY)] = tr (v exp (i§^A^j ^ . For contractible loops 7.;, 

we expect (W[yi] W[j 2 ] ■ ■ ■ W[7n]) = d n . At g = 0, this is 
obtained with d = 2 by direct calculation. This is somewhat 
surprising since one might have naively anticipated power-law 
correlation functions for fields governed by the gapless action 
(|9}. However, a peculiar feature of this non-relativistic action 
is that many equal-time correlation functions are short-ranged 
in space. For instance, dropping cubic and quartic terms, 
(B Q (x,0)B b (0,0)) - g 2 V 2 5W{x)5 ab . On the other-hand, 
we do expect power-laws to show up in non-local correlation 
functions such as: 

(tr (E t {x) V^ A 25,(0) Ve^ A )) ~ ±- % (12) 

where 7y 2 is an 0(n) loop model exponent defined after (|6j. 

Discussion. As we showed in the first part of this paper, 
the family of Hamiltonians Q is a fixed line of critical states 
which are most naturally described in the language of fluctu- 
ating unoriented loops. The Hamiltonian is defined by topo- 
logical relations which the loops must obey; these relations 
can be viewed as obtained by relaxing one of the relations 
(Jones-Wenzl) related to doubled SU(2)k Chern-Simons the- 
ory, thereby leaving only d-isotopy. In the second half of the 



paper, we found another fixed line of critical states, the non- 
relativistic SU(2) gauge theory fl5J, whose Wilson loop rep- 
resentation is also in terms of fluctuating unoriented loops. 
These are the two main results of this paper. It seems improb- 
able that two similar remarkable occurrences could be unre- 
lated, so we conjecture that the gauge theory (|9j controls the 
infrared behavior of the universality class containing the mi- 
croscopic model Q. 

There are a number of obvious generalizations and open 
questions which we hope to address later, such as the instabil- 
ities of these critical theories, especially at the magic values 
d = 2cos(7r/(fc + 2)); the addition of Chern-Simons terms 
to the action; correlation functions of and the precise rela- 
tionship between g and d; other gauge groups; complex d and 
0-terms in the action; the effects of instantons in (|9j; and the 
implications of our results for topological quantum computa- 
tion |d[ll. 
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